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ABSTRACT 

Aims. We discuss the importance of pure hydrogen white dwarf atmosphere models with Ly-a far red wing opacity in the analysis of 
the white dwarf cooling sequence of the globular cluster NGC 6397. 

Methods. Our recently improved atmosphere models account for the previously missing opacity from the Ly-a hydrogen line broad- 
ened by collisions of the absorbing hydrogen atoms with molecular and atomic hydrogen. These models are the first that well repro- 
duce the UV colors and spectral energy distributions of cool white dwarfs with r cff < 6000 K observed in the Galactic Disk. 
Results. Fitting the observed F814W magnitude and F606W - F814W color we obtained a value for the true distance modulus, 
fi = 12.00 ± 0.02, that is in agreement with recent analyses. We show that the stars at the end of the cooling sequence appear to be 
~ 160 K cooler when models that account for Ly-a opacity are used. This indicates that the age of NGC 6397 derived from the white 
dwarf cooling sequence using atmosphere models that do not include the correct Ly-a opacity is underestimated by ~ 0.5 Gyr. 
Conclusions. Our analysis shows that it is essential to use white dwarf atmosphere models with Ly-a opacity for precise dating of old 
stellar populations from white dwarf cooling sequences. 
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1. Introduction 

Globular clusters are the oldest known stellar systems in our 
Galaxy. They contain a large population of very cool and very 
old white dwarfs, remnants of stars that formed and died through 
the history of a cluster. Such populations have been investi- 
gated in the two nearest globular clust e rs: M 4 dHansen et al.L 
200 2t iRichter et aUl2004t Ha nsen et all I2004T) and NGC 6397 
dHansen et all l2007t iRicher et all I2006I) . Recent data on NGC 
6397 show a clear white dwarf cooling sequence and a low lumi- 
nosity cut off rel ated to the finite age o f the cluster. On the basis 
of this sequence, lHansen et al] (2007) derived the true distance 
modulus fio = 12.02 + 0.06, the reddening toward NGC 6397 
E(F606 W - F8 14) = 0.20 + 0.03, and its age T c = 1 1 .47 + 0.47 
Gyr. The detailed and extended analysis presented in that publi- 
cation gives the first precise and independent dating of a globular 
cluster based on its white dwarf cooling sequence. 

Reliable atmosphere model s need to be used in such analyses 
dFontaine et al.L 120011: iHansenl 1 19991). Models current ly in use 



eg. Bergeron et al. (2001); Sau mon & Jacobso: 



i currently 



Hansen 

119981) : IBergeron et al.l d 19951) ), although valuable in the anal- 
ysis of cool white dwarfs (eg. Berg eron et alj d2001l) ). have a 
well known shortcoming. They cannot correctly reproduce the 
UV flux ( A < 5000 A) of white dwarfs with r eff < 6000 K 
dBergerorl 1200 U IBergeron et ail l200ll 119971) . This problem 
has been solved recently by introducing the far red wing pro- 
file of hydrogen Ly-a into the modeling. The line profile is 
formed by collisions between absorbing hydrogen atoms and 
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the atomic and molecular hydrogen, and when introduced into 
the modeling it allows for an accurate reproduction of the spec- 
tral energy distributio n of co ol white dwarfs at all wavelengths 
dKowalski & Saumonl 120061) . The synthetic spectra resulting 
from our improved pure hydrogen atmosphere models deviate 
significantly from the previous ones for wavelengths shorter than 
1 //m (Fig . |TJ . This indicates that the Ly-a opacity also affects the 
flux in filters F606W and F8 14W used to obtain the photometric 
data on the white dwarf cooling sequence of NGC 6397. 

In this paper we discuss the importance of our new white 
dwarf atmosphere models with Ly-a opacity in the analysis of 
white dwarf cooling sequences. Fitting the theoretical to ob- 
served white dwarf cooling sequence of NGC 6397 we derived 
the true distance modulus to the cluster. We show that its value 
is consistent with that obtained with old models. We also predict 
that the coolest white dwarfs observed in NGC 6397 are cooler, 
when the atmosphere models with Ly-a opacity are used in the 
analysis. This implies that t he cluster age derive d from the white 
dwarf cooling sequence by lHansen et alj d2007l) is probably too 
low. Comparing the temperature of the white dwarfs at the cut off 
of the white dwarf cooling sequence in NGC 6397 and the galac- 
tic disk we show that the local population of disk white dwarfs 
is much younger than that observed in the globular cluster. 



2. Cool white dwarf atmosphere models 

lHansen et all d2007l) and lRicher et alj d2006l) in their analysis of 
the white dwarf coo ling sequence o f NGC 6397 used the at- 
mosphere models of Berger on et alj d2QQl). However, there is 
a well known shortcoming of these models, as they overestimate 
the flux at short wavele ngths (A < 5000 A, see Ber geron et al.l 
d 1997 ): Bergeron (2001)). We have corrected this deficiency by 
introducing the opacity of the far red wing of the Ly-a broad- 
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Fig. 1. Synthetic spectra of pure hydrogen atmosphere models 
with r eff = 3000 K, 4000 K, and 5000 K (from top to bottom), 
and fixed gravity log g — 8 (cgs) calculated with (solid) and with- 
out (dotted) L y-a far red wing opaci ty. Th e dotted lines mimic 
t he m odels of Bergeron et a lT d 1 995b and [Saumon & Jac obsonl 
(119991) used bv lHansen et aU(l2007t) and iRicher etail (120061) in 
their analysis of the white dwarf cooling sequence of NGC 6397. 



ened by collisions of the absorbing hydrogen atoms w ith atomic 
and molecular hydrogen (Kowalski & Saumon, 2006). Our new 
models allow for fits to the entire spectral energy distribution 
of white dwarfs observed in our Galaxy, including the fit to 
the near-UV HST spectrum of the white dwarf star BMP 4729 
dKowalski & Saumonll2006l:lKoester & Wolfltl2000l) . and the re- 
production of the observed white dwarf co lor sequences below 
~ 6000 K (iKowalski & Saumoni 120061) . 
In Fig. Q] we present the impact of the Ly-a opacity on the 
synthetic spectra of cool white dwarfs. The flux suppressed by 
the Ly-a opacity at short wavelengths is redistributed towards 
longer wavelengths (well beyond A ~ 1 yum). This causes a non 
negligible increase of the flux at visible and near infrared wave- 
lengths. As a result, the synthetic photometry in the F606W and 
F814W filters is also affected. This is shown in figure [2] where 
we present the color-magnitude diagram for a cooling sequence 
of a white dwarf of mass 0.5 Mq computed using models with 
and without the Ly-a opacity. Both curves deviate significantly 
from each other for T a s < 5000 K. As Ly-a opacity is stronger 
at short wavelengths, the effect is stronger in the F6Q6W filter 
(larger horizontal shift in Fig. 2). The new models predict a turn 
towa rds the blue, caused by CIA opacity from molecular hydro- 
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Fig. 2. Color-magnitude diagram of the white dwarf cooling se- 
quence. The mass of the star is 0.5 Mq, and the effective temper- 
ature runs from 6000 K (top) to 3000 K (bottom). Crosses indi- 
cate the T e ff with 500 K decrements along the curves. Predictions 
of both models with (solid) and without (dotted) Ly-a opacity 
are shown. On the vertical axis the absolute magnitude is indi- 
cated. 
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gen (Sau mon & Jacobsonl TT999; Hansen, 1998; Berg eron et al 



119951) . at lower Teg. This should have an impact on the astro- 
nomical parameters of NGC 6397 derived from the analysis of 
its white dwarf cooling sequence. 



3. Analysis of the data 

3. 1 . The fit to the white dwarf cooling sequence 

To investigate the importance of pure hydrogen white dwarf at- 
mosphere models with Ly-a- opacity in the analysis of the white 
dwarf cooling sequence of NGC 6397 we performed a least 
square 2 D fit in F814W and F606W-F814W to th e photomet- 
ric data of lHansen etail (l2007l) : lRicher et all (I2006I) : their Fig. 6 
and 4 respectively. The fitting parameters were the reddening and 
the true distance modulus. The fit was performed using two grids 
of atmosphere models: those with Ly-a opacity and those with- 
out. The second set mimics the models used in the analysis of the 
cluster whi te dwarf cooling se quence conducted bv lHansen et alj 
(2007) and lRicher et al. (2006). Following those authors, in our 
derivation of the theoretical cooling sequences we assumed the 
mass of the white dwarf to be 0.5 Mq. Our fits to the data are 
given in Fig. [3] We obtained perfect fits with both the o retica l 
cooli ng sequences, as i s also the case in lHansen et al.l ([2007) 
and iRicher et alj (I2006I) . In addition to that work we show, for 
the first time, the T s g of cool white dwarfs in the globular clus- 
ter. Because of the difference in localization of the turn off visi- 
ble in both theoretical sequences (Fig. [2]), both fits give different 
estimations of T e g for the c ool end of the obse rved sequence 
localized at F814W- 27.6 (lHansen et al.L I2007I) . The effective 
temperature estimated on the basis of atmosphere models that 
account for Ly-a opacity is significantly smaller, indicating that 
the stars at the cut off of the observed white dwarf cooling se- 
quence are cooler and older. We will return to this problem in 
the next section. 
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Table 1. Effective temperature and ages of white dwarfs at the 
end of the cooling sequence of NGC 6397. 

mass(M ) r eff (K) age (Gyr) T eS (K) age (Gyr) 



F606W-F814W 



Fig. 3. Color-magnitude diagram of the white dwarf cooling se- 
quence observed in NGC 6397. The lines represent the fits to 
data using modeled cooling sequences for a white dwarf of 
mass 0.5 Mq. The effective temperature runs from 6000 K (top) 
to 3000 K (bottom). Small crosses indicate the r e jf with 500 K 
decrements along the curves. Models with (solid) and without 
(dotted) Ly-a opacity are shown. 



3.2. Reddening and the true distance modulus 

As indicated in Figure [2] there is a significant variation in the 
modeled maximum value of the F6Q6W - F814W color for a 
white dwarf of a given mass. In the present case, both model 
curves deviates by ~ 0.1 magnitude, which has a direct impact 
on the derived value for the reddening towards NGC 6397. With 
models that do not account for the Ly-a opacity we obtained 
E(F6Q6 W-F8 UW) = 0.1 8 +0.2 with a goo d agreement with the 
value 0.16 obtained bv lHansen et"aT1(l2007l) . Models with the Ly- 
a give a lower value of 0.12 + 0.02. This result, when converted 
from the ACS photomet ric system to the g round-based colors 
using the prescription of Siri anni et al.l d2005l) gives E(B - V) = 
0.1 3 + 0.03. This value is lower than E(B - V) = 0.18 obtained 
by iGratton et al.l d2003l) using a fit to the Main Sequence stars. 
However, the reddening obtained by fitting the M = 0.5Mq the- 
oretical cooling sequence to the photometric data is expected to 
be underestimated. This is because of the variation in the average 
mass of a white dwarf a t the end of the coo ling sequence. This 
is discussed in detail by lHansen et al.l d2007l) . The final redden- 
ing derived in that paper is E(F6Q6W - F814W) = 0.20 + 0.3. 
This is ~ 0.04 dex larger than the value they obtained using the 
M = 0.5Mq theoretical cooling sequence. Applying such a shift 
to our result, we obtain E(F6Q6W - FUAW) = 0.16 + 0.02 and 
E(B - V) = 0.17 ± 0.03 respectively. This va lue is in good agree- 
ment with the result of lGratton et all (120031) . 

For the true distance modulus, in both fits we obtained a con- 
sistent val ue of 12.00 + 0.02 , which agrees with the value ob- 
tained byEaniinitin (|2007i), ,u = 12.02 + 0.06. Our estimate of 
the true distance modulus was obtained assuming a single white 
dwarf mass of M — 0.5Mq and t he relation between redden- 
ing in both filters Ash = 0.65A606 dSirianni et al.L 120051) . When 
we use the M = Q.6Mq mass star instead, the resulting distance 
modulus is lower, 1 1.85 + 0.02. 



models with Ly-a 
0.5 3180 ± 80 9.8 + 0.4 

0.6 3210 + 80 11.1 + 0.3 



models without Ly-a 
3340 ± 80 9.2 ± 0.3 
3380 ± 90 10.6 ± 0.3 



4. r e ff at the edge of the white dwarf cooling 
sequence and the age of NGC 6397 

The proper derivation of the age of NGC 6397 from its white 
dwarf cooling sequence is a difficult task and is beyond the scope 
of this paper. However, as use of the new atmosphere models 
shifts the effective temperature of the coolest white dwarfs by a 
few hundred degrees (Fig. [2]), it is important to show how large a 
bias this effect would produce on the derived age of NGC 6397. 
To check this we cond ucted a simple exercise, similar to that 
of lHansen et alj d2007f) . appendix A. Following the authors we 
assumed that the observed cut off of the white dwarf cooling se- 
quence takes place at F814W = 27.6 + 1. Then assuming the 
white dwarf masses to be 0.5 Mq and 0.6 Mq respectively we 
deriv ed the r e ff and, usin g t he same white dwar f cooling mod- 
els of iRicher etaTI d2000l) as lHansen etal.1 d2007l) , we derive the 
cooling times of the white dwarfs at the cooling sequence cut 
off. The results are summarized in TableQ] The T e ff and the ages 
derived using atmosphere models with and without Ly-a opacity 
differ by ~ 160 K and ~ 0.5 Gyr respectively, and these differ- 
ences are insensitive to the assumed white dwarf m ass. This sug- 
gests t hat the age of the globular cluster derived in lHansen et al.l 
(2007) can be ~ 0.5 Gyr too low because of the white dwarf at- 
mosphere models that do not account for Ly-a red wing opacity. 

We also conducted the same simple analysis for the white 
dwarf cooling sequence observed in the Milky Way o b tained 
from three samples of white dwar fs dBergeron et al.L 2001; 
iKilic et ail 120061: lHarris et"a"Q |2006[). To do th i s we u sed pub- 
lished color sequences of Kowalski & Saumon (2006). The ob- 
tained value of colors and T e ff at the cut off are given in Table 
[2] Our analysis shows that the coolest white dwarfs observed in 
the Galactic Disk have r e jf ~ 3800 K. This indicates that the 
coolest stars in the local population of white dwarfs are warmer 
by ~ 600 K than the coolest white dwarfs obser ved in NGC 
6397. According to white dwarf cooling models of IRicher et al.l 
(2000) such a differ ence in effective temperature results in ~ 2 
Gyr difference between the age of the white dwarfs in NGC 6397 
and the local whi te dwarf population ( Fig.[4]i. This is consistent 
with the result of Fontai ne et al.l d2001 ), who, based on the anal- 
ysis of the white dwarf luminosity functions, concluded that the 
Galactic Halo is 2 - 3 Gyr older than the Galactic Disk. 



Table 2. Effective temperature of the white dwarfs at the end of 
the cooling sequence observed in the Galactic Disk. 



color 


value 


T* (K) 


B-V 


1.4 + 0.1 


3750 ± 250 


V-I 


1.4 + 0.1 


3800 ± 300 


u-g 


2.6 ± 0.1 


3750 ± 150 


g-z 


1.9 + 0.1 


3850 ± 150 
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Fig. 4. White dwarf cooling times vs. T e g for a white dwarf of 
mass 0.5 Mp) (soli d line) and 0.6 Mq (dotted line). Models of 
iRicher et al.l(l2000l) . 



5. Conclusions 

Observation of the white dwarf cooling sequences in the glob- 
ular clusters opens the possibility to date them more precisely. 
However, there are many factors that affect the models applied 
to white dwarf cosmochronology. Here we have discussed the 
impact of white dwarf atmosphere models on the derivation of 
the basic astronomical properties of NGC 6397 from the white 
dwarf cooling sequence itself. Our most recent set of pure hydro- 
gen models, that corrects for the flux discrepancy at short wave- 
lengths, results in lower reddening toward the cluster, and older 
derived ages. By deriving the effective temperature of white 
dwarfs at the end of the white dwarf cooling sequences of NGC 
6397 and the Galactic Disk we have also shown that the coolest 
white dwarfs observed in the globular cluster are ~ 600 K cooler, 
and therefore ~ 2 Gyr older. More sophisticated analysis that in- 
cludes the construction of the white dwarf luminosity function 
is necessary for correct derivation of the cluster age from the 
white dwarf cooling sequence; however, our analysis shows that 
atmosphere models with Ly-a opacity are necessary to reliably 
estimate the ages of old stellar populations from the white dwarf 
cooling sequences Q. 
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